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Abstract 



In minimal SUSY GUT models the .R-parity breaking terms are severely constrained 
by SU(5) gauge invariance. We consider the particular case where the explicit -R-parity 
breaking occurs only via dimension 2 terms of the superpotential. This model possesses 
only three R-parity breaking parameters. We have studied the predictions of this model 
for the neutrino masses and mixing angles at the one-loop level within the framework of 
a radiatively broken unified supergravity model. We find that this model naturally yields 
masses and mixing angles that can explain the solar and atmospheric neutrino problems. 
In addition, there are regions in parameter space where the solution to the solar neutrino 
puzzle is compatible with either the LSND result or the existence of significant hot dark 
matter neutrinos. 



Chapter 1 
Introduction 



The standard model of elementary particle physics (SM) is in very good agreement with 
all presently available data. Nonetheless, it suffers from various theoretical shortcomings, 
the most severe of which is the hierarchy problem]!]]. 

In supersymmetric models the cancellation of quadratic divergences is guaranteed and 
thus any mass scale is stable under radiative corrections. Supersymmetry (SUSY) implies 
that any fermion (boson) is accompanied by a bosonic (fermionic) superpartner with 
the same mass and transformation properties under the gauge symmetry 0. The most 
economical candidate for a realistic model is the minimal supersymmetric extension of the 
SM (MSSM). In addition to the superpartners for all SM particles it contains two Higgs 
bosons (required to give mass to up and down-type fermions and to cancel the triangle 
anomalies arising form the fermionic partners of the Higgs bosons) but no other particle. 

The most general superpotential invariant under the gauge symmetry can be written as 

W = ly^LrLjEi + yf^UQ^l-ylHQ^ (1.1) 

where the supermultiplets are denoted by a hat. The left-handed lepton supermultiplets 
are denoted by Lj (i — 1,2,3) and the Higgs superfield coupling to the down-type quarks 
is denoted by L . Throughout this paper, we use the notation i, j, k = 1,2,3 and I, J, K = 
0, 1,2,3 and we sum over twice occurring indices. Note that LiLj = e a bLfL b j = —LjLi 
(a, b = 1, 2, 3 are the SU^)^ indices) and thus yf Jk = —yj Ik , 

Let us first determine the meaning of the various terms of eq. |lTTi Here, y^- fc , y^ k and 
y^ k denote the lepton, down-type and up-type Yukawa couplings, respectively, and fio is 
the Higgs mass parameter. However, in contrast to the SM the MSSM allows for baryon 
[lepton] number violating interactions y? k [y^ k and /ij]. 

These couplings are constrained from above by experiment. The most model independent 
constraints can be obtained from collider experiments [§] or neutrino physics [@,||]. It 
turns out that the individual lepton and baryon number violating couplings only have to 
be smaller than 0(1CT 3 ~ few x 1CT 1 ). Thus, the i?-parity violating couplings need not 
be more suppressed than the lepton and baryon number preserving Yukawa couplings. 
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The exception is the constraint on the coupling 

y?2i & io~ 7 (1.2) 

from heavy nuclei decay f|, but even this constraint becomes somewhat less impressive 
if compared to y^ u ~ 3 x 10~ 5 . Somewhat stronger constraints can be derived from 
cosmology 0. 

Thus, it may be premature to conclude from our negative experimental search that im- 
parity is a good approximation or even exact symmetry of nature. 

However, the experimental exclusion area can be strongly enhanced by imposing theo- 
retical constraints. In the minimal SU(5) SUSY-GUT model, the right-handed leptons, 
the right-handed up-type quarks and the left-handed quarks are embedded in a ten di- 
mensional representation, 10j = Ef © Uf © Qi. The right-handed down-type quarks and 
the left-handed leptons are embedded in a five dimensional representation, 5j = D i © Lj. 
The two Higgs doublets are embedded together with two proton decay mediating colored 
triplets, T and D , in five dimensional representations, 5 = D Q © L Q and 5 = T © H . In 
this model, both the lepton and the baryon number violating interactions arise from the 
term 

Wgut = §2/ijfc5j5jl0jfc , (1.3) 

where the boundary conditions at Mqut are given by y^ k = yf k ^ = y® k = y-ijk^ Thus, in 
general the baryon and lepton number violating couplings are correlated in SUSY GUT 
models. This leads to very strong constraints on any from proton stability which are 
much stronger then the constraint in eq. Q|8fl. (Note, that it does not constraint the 
coefficients of dimension 2 terms, 

These much stronger constraints in the framework of SUSY GUT models are the reason 
why in the MSSM any lepton and baryon number violating interaction is eliminated by 
imposing a discrete, multiplicative symmetry called i2-parity[[| 

R p = (-lY s +* B +\ (1.4) 

where S, B and L are the spin, baryon and lepton numbers, respectively. Aside from 
the long proton life-time, i?-parity conserving models have the very attractive feature 
that the lightest supersymmetric particle (LSP) is stable and a good cold dark matter 



candidate |T0[. However, while the existence of a dark matter candidate is a very desirable 
prediction, it does not prove i?-parity conservation and one should keep an open mind for 
more general models. 

In this paper, we will investigate the scenario where i?-parity is broken explicitly via 



Hi 7^ 0|1. In particular, we compute the predictions for the neutrino masses and mixing 



1 It is clear that the predictions for the down-type quark masses of the first two generation in the 
minimal model are off by factors of 0(3). Thus, any more realistic model has to be more complicated. 
While it is easy to reconcile the mass predictions with experiment, e.g. by introducing higher dimensional 
representations or higher dimensional operators, it would be very hard to explain why either the lepton 
number violating Yukawa couplings or the baryon number violating couplings are suppressed by many 
orders of magnitude. 
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angles in the frame-work of a SUSY-GUT model with radiative electroweak symmetry 
breaking. 

Our paper is organized as follows: in section 2 we present the model at tree-level and our 
notation and conventions. In section 3 we present our numerical analysis of the one-loop 
radiatively corrected neutrino/neutralino mass and mixing matrix. Our conclusions are 
presented in section 4 and our results for the renormalization group equations (RGEs), 
the expressions for masses and vertices, and the formulae for our one-loop results are 
relegated to four appendices. 
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Chapter 2 

The Model at Tree-Level 



In the introduction we have argued that explicit i?-parity breaking Yukawa couplings are 
strongly constrained by the long proton life-time if the model originates from a SUSY- 
GUT model. One alternative is to break i?-parity spontaneously at the electro-weak scale. 
This possibility has been studied extensively [[|5|,[p^],|13|1. The main feature of these mod- 



els is the existence of a Goldstone-boson, the Majoron, associated with the spontaneous 
breaking of a continuous symmetry, Lepton number conservation. The simplest model 



where i?-parity is broken spontaneously by a sneutrino VEV[T2[] is phenomenologically 
ruled out. This leads to the introduction of several SM singlets and new Yukawa- type 
couplings |jl3|. Another possibility is the explicit -R-parity breaking via soft terms of the 



superpotential only. These terms are unconstrained by the proton life-time and can arise 
naturally in SUSY-GUTs. 

Let us consider, for example, a SUSY GUT model based on SO(10). Here, the Higgs fields 
are embedded in a 10 = 5 © 5o while the down-type fermions are embedded together 
with all other fermions of one generation including the right-handed neutrino in a 16 = 
10 © 5 © 1. The only renormalizable Yukawa couplings in SO (10) can be written as 
Wqut = VijlO © 16j © 16j. This automatically implies the absence of the couplings y^ k1 



y^ k and y^ k . The only non-zero components are from now on abbreviated by y^ = y^ 
and yfj = y^. 



Furthermore, the Higgs/slepton mixing terms vanish ( i.e. fii — 0). However, in general 
this term can be generated dynamically via W m %x = 2/ilO © 16j © 16h if SO (10) is broken 
spontaneously to SU(5) by the right-handed neutrino-like component of the 16 dimensional 
Higgs field (16//) ^ 0. Why the resulting mass terms ^ = Aj(16#) and fi are of the 
order of m z rather than the scale of gauge unification, Mqut — 10 16 GeV is not clear 
but possible mechanisms are known |T4|. It is not implausible that the couplings yi arise 
from non-renormalizable terms and are 0{m^/2/M-p) (here, m.3/2 — fewxlOO GeV is the 
gravitino mass and Mp ~ 10 19 GeV is the Planck mass). Whatever the solution for the 
\i problem, it is expected to affect the fii in an analogous fashion. Thus we can expect 
naturally that /ij = 0(/io)[] with a possible suppression factor of (16#)/Mp ~ 10~ 3 . In 

1 the only plausible alternative would be [ii = Mqut which is phenomenologically ruled out. 
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this model, the analogous baryon number violating Higgs triplet / down- type quark mixing 
is suppressed automatically if D (together with T) acquires a mass O(Mgut)- Such a 
large mass for the Higgs triplets is required in any SUSY-GUT model with or without 
i?-parity. It is a very severe problem[|15[] and has no completely satisfying solution to 
date. However, once this problem, common to all SUSY-GUT models, is solved there are 
no additional baryon number violating terms in the low energy theory even if we include 
non-zero values of m. 

On the other hand, lepton number violation via Higgs/slepton mixing is potentially un- 
suppressed. This gives rise to neutrino masses and mixing angles and may provide an 
economical explanation for the solar neutrino deficiency ||16|| and to solve the atmospheric 
neutrino puzzle Hl7j . These possibilities are well within present constraints and deserve to 
be studied. Thus, we will focus our attention to the scenario where = but ^ ^ 0. 



2.1 The SUSY-GUT framework 

In any realistic model SUSY has to be broken at the electro-weak scale given by the 
mass of the Z boson, m z = 91.187[|T^] in order to explain why no superpartner has been 
detected to date. Here, this is done by including explicit soft SUSY breaking terms in the 
potential 

V soit = \A\^UE) + A%l Q QiD) - AVHQiU? - BfLjH + H.c." 

+ E MJL^A + E M aXaXa + H.C. , (2.1) 

m,n a 

where <p n = Lj, Ei,Qi,Ui, Di, H stands for all the sfermion fields and the Higgs doublets. 
The gauginos are denoted by A a [the indices a,b = 1,2,3 refer to the U(l), SU(2)^ and 
SU(3) C gauge symmetries, respectively; the gauge indices are suppressed]. 

The large number of free soft SUSY breaking parameters can be reduced significantly by 
making certain assumptions about their origin. In minimal N — 1 supergravity (SUGRA) 
models the breaking of SUSY is thought to occur in a "hidden" sector i.e. a sector of 
the theory that couples to the standard model particles only via gravity. The soft SUSY 
breaking terms arising from gravitational coupling to the "hidden sector" are assumed to 
be universal at the Planck scale Mp ~ 10 19 GeV. As a result, there are only four indepen- 
dent parameters. These are the A (B ) parameter multiplying the tri-linear (bilinear) 
part of W, the universal scalar mass m and the universal gaugino mass parameter m^, 
i.e. 

Auk = A y IJk , Af jk = A yf jk , A tj = A y^. , 
Bj = B L jii j = B fMi , 

M a = m 1/2 . (2.2) 

By assuming universal soft SUSY breaking parameters at Mp we have drastically reduced 
the number of free SUSY breaking parameters of the theory. This universality of the soft 
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SUSY breaking parameters is broken at the electro-weak scale due to a renormalization 
group evolution of the individual parameters^ (the (3 functions for the general model with 
i?-parity breaking are listed in appendix A). It is one of the great successes of this scenario 
that the low energy value of m\ obtained from renormalization group evolution is indeed 
negative, giving rise to spontaneous electroweak symmetry breaking while all other scalar 
mass parameters remain positive over a large region of parameter space |E0 . 



2.2 Minimizing the Higgs Potential 

Before we can investigate any other sector of the model we first have to minimize the 
Higgs potential. The fundamental difference of our model to the MSSM with conserved 
i?-parity is the fact that the Higgs fields mix with slepton fields. Thus, we effectively have 
to minimize a five Higgs doublet model. In particular, all the mass eigenstates and the 
vacuum expectation values will be a linear combination of Higgs fields and slepton fields. 
The tree-level potential can be written as 

V = (ti 2 + m%)rfH+(wj + mljL\L J -B?(L I H + tt.c) 

+ 9 -^- (h*H - IP/) 2 + £ \H%\ 2 . (2.3) 

The soft SUSY breaking parameters are obtained by numerically solving the renormal- 
ization group equations (RGE) given in Appendix A with universal boundary condi- 
tions [eq. |2~2]1 . For a qualitative understanding we present the result for tan (3 = 20 and 
m t = 175 GeV. 



Bjj = [B + 0.25m 1/2 - 0.354) S u + AB LJ , 
™L = (m 2 + 0.5mj /2 ) 5 U + Am 2 Lij , (2.4) 

The non-universal terms are generated through the down-type Yukawa couplings. For a 



qualitative understanding it is convenient to neglect these terms|2l|. They are to a good 
approximation independent of m t assuming the top Yukawa coupling is not too close to 
the Landau pole and have a trivial tan (3 dependence 

10 -4 

AB m ~ — — (2.2m 1/2 - 1.24,) , 

COS ID 

AB l3 ~ ^^(0.127^/2-0.34), 
Am 2 Lm ~ ( 3 - 7m o + n.8ml /2 + 1.0^ - 3.6A m 1/2 ) , 

COS l) 

Am!,, - -y^yf k (0.9ml + 0.25m 2 1/2 + 0.28A 2 -0.21A m 1/2 ) , (2.5) 



2c 



'Strictly speaking, this scenario assumes the absence of any intermediate scale between Mp and m z 
since significant non-universal terms can arise already in the relatively small range between Mp and 
Mq\jt due to large group coefficients of a unified groupjl^. However, these effects are strongly model 
dependent and shall be neglected here. 
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with all other off-diagonal elements vanishing. We introduce the following notation for 
the individual components 

H + \ r _/(£!- «BJ)/V2 



The values of v and Vj are obtained by minimizing the potential in eq. p.3| numerically. For 
small i?-parity violating parameters we can also obtain very reliable analytic expressions 
in the basis where y£ is diagonal by expanding in powers of /i«//io 

• q/3 2i? 00 /i 

Sm213 = m L+ < + 2,g - (2 ' 7 > 

W/3 = (2 . 8 ) 

Vi_ B (ii) tan/3 - /i 

+ ^ cos 2/3 

with the convention that indices in braces are not summed over. In order to stay as close 
to the notation of the MSSM as possible we define 

(H°) (L°) 

v = — , vi = — j=r , v = y/vivi , and tan (3 = v/v, (2-10) 
v2 v2 

and we parameterize the VEVs in terms of spherical coordinates 

tan 04 = — , tan0' = — ^— , tantf = — ^— . (2.11) 

t>2 V\ COS 6>3 t>o COS 02 

Analogously, it is convenient to parameterize the .R-parity breaking mass parameters in 
terms of three mixing angles 

tanfl 3 = — , tan# 2 = — — , tan^ = — ^ . (2.12) 

[12 yUl COS 6 3 yU cos 6 2 



and [i = JJltJii- The potential in eq. |273] is minimized by an iterative procedure using 
the analytic solution for tan#i = as our initial values. This procedure also works 
surprisingly well for tan#i > 1. However, in order to obtain qualitative understanding of 
the results it is instructive to investigate the potential analytically. Let us for the moment 
neglect the effects of the down-type Yukawa couplings on the running of the soft SUSY 
breaking parameters. Then the conditions m 2 Lij oc 5jj and Bj oc /i/ would continue hold 
even after renormalization group evolution down to the electro-weak scale. Let us make 
a rotation on the Higgs/lepton superfields Lj — > U\jLj with the unitary SU(4) matrix 
defined by: 



pnUu = V// 2 (l, 0, 0, 0) . (2.13) 

In this basis there is no .R-parity violation in the Higgs potential of the model. Instead 
there are -R-parity violating Yukawa couplings. For a qualitative understanding this basis 
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is much more convenient ||1 1|| . In particular, we find that universal soft SUSY breaking 
parameter at m z lead to alignment of vi and fii which means 0i = 6[. It is clear that by 
solving the RGEs including the bottom Yukawa coupling, y^ oc mfttan/3, we have 

This result shows no explicit tan/5 dependence (for fixed Yukawa couplings). However, 
eq. |2.14| is too naive as we will see by considering the limit of large tan/3. In this case 



Bqq vanishes [eq. [2.7|| but B^ does not. This leads to an enhancement of the ratio Vi/v 
[eq. [2.9|| . In general, we find that eq. pT7^pT9| become unreliable if ABjj > B 00 . In this 



case v j and fij are completely misaligned and our numerical minimization procedure must 
be applied. 

Let us now turn to the complete particle spectrum of the Higgs sector. There are five 
CP-even neutral scalars, five CP-odd neutral scalars and eight charged scalars 

H° x = Ug{H°,L°) y (x,y = l,..,5), 
A° x = U^{A\B«) y (x,y = l,..,5), 

H± = UjtiH^Lf,^ (a,6=l,..,8). (2.15) 

This fields are mixed states of Higgs bosons and sleptons. The unitary matrices 
((f) = H°, A , H ) are obtained by diagonalizing the corresponding mass matrices given 
in Appendix B. We define the unitary matrices such that < m^. for i < j. Note, 
that the CP-odd neutral fields and the charged fields contain an unphysical mass-less 
goldstone-bosons, A\ and Hf, which are absorbed via the Higgs mechanism. These fields 
cannot occur as external fields but they do contribute in loops. 
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Chapter 3 

Neutrino/Neutralino Sector 



The i?-parity violating parameters /ij are expected to have the most noticeable effect 
on the neutrino/neutralino sector. In models with broken i?-parity the neutrinos and 
neutralinos mix via the mass matrix in eq. |B.2[ However, only one of the neutrinos, 



presumable u T , acquires a mass at tree-level contrary to the Majoron models where also 
m v is generated at tree-level^. Furthermore, the neutrino/neutralino mixing turns out 
to be rather small in the phenomenologically allowed region. Thus, it is natural to keep 
the conventional MSSM terminology and to refer to the three lightest (four heaviest) mass 
eigenstates as the neutrinos (neutralinos). 

In order to obtain the leading contribution to the masses of the two lightest neutrinos we 
have to compute one-loop radiative corrections to the mass matrix. 




= 1 



[l 2^ n + r e1]( p 2 )-i [l < n + r n *J(p 



Figure 3.1: The Lorentz invariant decomposition of the neutrino/neutralino self energy. 



The complete radiative corrections to the neutralino mass matrix in the MSSM with un- 



broken i?-parity has been presented in ref. [^2J and the leading order term for the neutrino 
masses has already been derived in ref. The one- loop corrected neutrino/neutralino 
mass matrix is given by 



n L(™ 2 ) - "m s L(V 



A=0 



(3.1) 



Here, the divergences are regularized by dimensional reduction|23|] where /ir is the renor- 
malization scale and d is the number of space-time dimensions. The superscript DR 
denotes the renormalized mass matrix in a minimal subtraction scheme obtained by set- 
ting A = 2/(4 — d) — 7£ + ln(47r) = where 7e is the Euler constant. The Lorentz 
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invariant amplitudes are defined in fig. |3.1| and the result for the one-loop self energies 
are presented in Appendix D[| 



For the mass eigenvalues it is usually sufficient to consider only the radiative corrections 
to the diagonal matrix elements of the propagator (i.e. n = m). However, in the case 
of mass-degenerate states at tree-level such as the two lightest neutrinos one finds the 
mass eigenvalues by rediagonalizing the one-loop corrected neutralino mass matrix. In 
addition, we are interested in neutrino oscillations determined by the mixing angles. Since 
we are only interested in the leading effects it would be sufficient to consider only the one- 
loop corrections to the v^-v^ mixing which is undetermined at tree-level. However, for 
a systematic treatment we evaluate the one-loop corrections to all matrix elements and 
rediagonalize the full 7x7 mass matrix. 

We define the one-loop mixing matrix, AZ such that the matrix m x o = (1 + AZ)J\A x ° lc (l + 
AZ^) is diagonal with the the one-loop radiatively corrected neutrino/neutralino masses 
m x o as the diagonal elements. The choice for the momentum at which the off-diagonal 
elements have to be evaluated is ambiguous. However, the effects of this ambiguity is of 
higher order since the off-diagonal elements are only important if m x o — m x o < 0(a/4-7r) 
and we simply choose p 2 = m 2 = (m 2 + m 2 )/2 

An Xm 

In a complete physical scheme any explicit dependence on gauge parameter or renormal- 
ization scale should cancel. A complete analysis would require the renormalization of the 
parameters v , g, g', tan (3, fi, M and M' at one-loop. However, we are only interested in 
the leading contributions to the neutrino masses and mixings. Thus, the one-loop correc- 
tions are only relevant for the two lightest neutrino masses. We have checked that the 
radiatively generated masses of the lightest two neutrinos are indeed gauge independent 
and renormalization scale independent. For the third neutrino (which acquires mass at 
tree-level) the scheme-dependence is of higher order and can be neglected. 



3.1 The neutrino mass spectrum: A general Scan 

In order to see what a typical neutrino mass spectrum looks like we first perform a general 
scan over the entire parameter space. We impose the standard assumption that all the 
soft SUSY breaking parameters are universal at Mgtjt- The spectrum then depends on 
four SUSY parameters, tan/3, A Q , mi/ 2 and fi (here, the value of mo is determined by 
imposing radiative electro- weak symmetry breaking and B is replaced in favor of tan /3) 
and in addition on three .R-parity violating angles tan#j (i = 1, 2, 3). 

From our Appendix C we see that the lepton current coupling to the W~ can be written 

1 It was shown in ref [pi| that the tadpole diagrams have to be included for a gauge-independent 
definition of the running masses. Here, we neglect the tadpole diagrams as they do not contribute to m„ e 
and m v . 
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as 



j; = ( L0 - + R °™) x * • (3 - 2) 

Here, \/20^ x /g is the analog to the CKM matrix of the quark sector. Note that in our 
notation the conjugate left-handed lepton doublets form the right-handed component of 
the charginos and vice versa. This unconventional arrangement was neccessary in order 
to have a unified notation for neutrinos and neutralinos and for charged leptons and 
charginos. 

The one-loop radiatively corrected interaction marices L and R are obtained from 
eq. |C14| by replacing Z — > (1 + AZ)Z. The analogous one-loop corrections to the charged 
lepton mixing matrix is suppressed by the inverse power of the lepton masses and can 
savely be neglected. 

There are two fundamental difference to the quark sector: (i) the existence of a left-handed 
and a right-handed current and (ii) the matrices \/20^ x /g (P = L, R) are not unitary. 
However, it turns out that only the case tandi <d 1 is phenomenologically interesting. 
In this case, the 3x3 submatrix <C 1 = 1,2,3) can be neglected. The 3x3 
submatrix Of, is to a very good approximation unitary and equivalent to the leptonic 
CKM matrix. Furthermore, in order to adopt a more standard notation we define 



V2 
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sm9 e ^ = — O* , i^j. (3.3) 



Strictly speaking, this is only correct in the case of two generation mixing or for small 
mixing angles, but the accuracy is still better than the experimental uncertainties and 
will be sufficient for our purposes. 



It is somewhat problematic to present the result of different diagrams given in eq. p. 2 



and p.l| because the mass eigenvalues of a sum of two matrices is not equal to sum of 
mass eigenvalues of two matrices. Thus we will only consider the full result obtained 
by summing over all one-loop diagrams. However, we would like to emphasize that the 
only significant diagrams are the ones envoi ving the down-type quarks/squarks and the 
charginos/charged Higgs bosons. It turns out that for small values of tan f3 the down-type 
quark/squark loops dominate due to the large Yukawa couplings and a color enhancement. 
For large tan/3 there are large sneutrino VEVs even in the basis where /xj oc (1,0,0,0). 
Here the dominant effects arise from the chargino/charged Higgs loop with sneutrino VEV 
insertions rather than from i?-parity violating Yukawa couplings. 



In fig. |3.2| and |3.3| we show the dependence of the three neutrino masses and the e-/x and 



jjtr-r mixing angles as a function of all seven parameters. Our master set of parameters 
is: A = 0, m-i/2 = 100 GeV, /x = 200 GeV, tan/3 = 20 and we assume maximum mixing, 
i.e. Hi = |/x(l, 1, 1, 1)/ or tan#i = \/3, tan# 2 = \/2, tan6* 3 = 1. Furthermore, we set 
m t = 175 GeV in all plots. 
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Figure 3.2: The one-loop radiatively corrected neutrino masses as a function of (a) tan#i, 
(b) tan $2, (c) tan 63 and (d) tan/5. Our master set of parameters is \x = 200 GeV, 
A = 0, mi/2 = 100 GeV and we chose tan 2 #i = 3, tan 2 # 2 = 2, tan 2 # 3 = 1, tan/3 = 20 
and m t — 175 GeV whenever these variables are not varied. 
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Figure 3.3: The one-loop radiatively corrected neutrino masses as a function of (a) mi/2, 



(b) Aq, (c) /i < and (d) > 0. Our master set of parameters is as in fig. K2 
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Figure 3.4: The one-loop radiatively corrected neutrino masses as a function of (a) the 
top quark mass and (b) the SUSY breaking scale parameterized by m^. We have set 
A = and /i = 2mi/ 2 . 

In fig. [37| (a)-(c) we see that 

fi = 1,2,3 for I = e , 

m Vl oc tan 2 Q, L , where < % = 1, 2 for £ = [i , (3.4) 

I % = 1 for £ = r , 



and 



sin 2 20 ei ^ oc tan 2 9 3 , 



sin 2 2^ T oc tan 2 6 2 . (3.5) 

as long as tan#j <ti 1. The dependence on tan/5 is more complicated. We see that there 
is a sharp rise in m Ur if the top Yukawa coupling is near a Landau pole ( i.e. tan/3 ~ 2). 
On the other hand, there is also an increase of all three neutrino masses with tan/3. Note 
that the explicit cos 2 (3 suppression of m Ur in eq. [B.3| is overcompensated by the tan 4 f3 
dependence of sin 2 (#i — 9[). 

In fig. |37^(d) no obvious proportionality can be established. We can see that a destructive 
interference of various one-loop diagrams can result in a vanishing m u for a particu- 
lar (large) value of tan/3. The same can happen to for some values of the SUSY 
parameters. 

We find that over most of the parameter space there is a hierarchy among the three 
neutrino masses. However, there are regions where the heaviest one-loop generated mass 
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can dominate over the tree-level mass. E.g. in fig. |3.3| (c) we see that m u and m Ur intersect 
for ji ~ —300 and —200 GeV. On the other hand, we find that the mass of the lightest 
neutrino m Ue <C m v , m Vr over the entire parameter space. 



In fig. |3.4j we present the masses and mixing angles as a function of the SUSY breaking 
scale parameterized by m 1( / 2 for A = and \i = 2m\/2- We find that all the masses are 
inversely proportional to the SUSY breaking scale (aside from singularities in the lightest 
mass eigenvalue due to some accidental cancellation among different diagrams). This 
decoupling is due to the fact that our model reduces to the SM in the limit of large SUSY 
breaking scales. 



Finally, in fig. we present scatter plots of 6 x 10 4 different sets of parameters. Again we 
assume maximum mixing and we have scanned over all SUSY breaking parameters in the 
range 2 < tan/5 < 40, \ A/m 1 / 2 \ < 3, 50 GeV < mi/ 2 , \fi\ < 500 GeV and 1 < \/j,/m 1 / 2 \ < 5. 
Typically, we find m v Jm Vr ^ 10 -4 and m Vf ,/m VT ^ 10~ 7 . The fact that the heavier one- 
loop radiatively generated mass and the tree-level mass are relatively close in magnitude 
can be understood from eqs. [2.14 and B.3 : 



m,. 



0i 



m,. 



0i - 6[ 



x O 



a. 



47T 



(3.6) 



[. Further- 



We see that the one-loop suppression is compensated by the smallness of Q\ 
more, we see that even in the case of maximum i?-parity violation the neutrino masses 
are below their experimental bounds|^5| over most of the parameter space. The reasons 



for this suppression was discussed in detail in ref. [pTT| 



A neutrino mass spectrum for different mixing angles can be obtained to a very good 
approximation by using the proportionality relations of eq. p.4| and |3~5. 



3.2 Solar and atmospheric Neutrino puzzles 

So far we have only considered the theoretically predicted neutrino spectrum. We have 
focused on the case of maximum mixing, from which the general case can be derived by 
rescaling the masses by appropriate powers of tan#j (i = 1, 2, 3). 

Before discussing any virtue of a given model we should assure that it satisfies all experi- 
mental bounds. There are various constraints on the SUSY spectrum. In particular, from 
LEP-1 we know that no superpartner (except a gaugino-like neutrino) can exist with mass 
below m z /2|1l8|. Stronger constraints on strongly interacting sparticle masses are obtained 



from experiments at FERMILAB||18||. All of these bounds can be evaded by scaling up the 
soft SUSY mass parameters just like in the models with unbroken i?-parity. Typically any 
model with m ,mi/ 2 ^ 100 GeV provides a sufficiently heavy sparticle spectrum. One 
can try to derive stronger constrains from virtual effects in processes as b — > 57, but it 
requires very strong assumptions on the squark mass matrices and possible cancellations 
among different contributions makes it hard to derive firm bounds that are much higher 



than those coming from direct particle searches p6 
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Figure 3.5: Scatter-plot of 6 x 10 4 models with maximal mixing in the (a) the m UT -m u ^ 
plane and (b) the m UT -m Ue plane. 
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In order to impose the experimental constraints on the Yukawa couplings on our model we 
have to rewrite the Lagrangian in a basis where the Higgs fields are the only ones with non- 
zero VEV. In this basis, there exist i?-parity violating Yukawa couplings proportional to 
the lepton masses y? k ~ yfk v i/ v o i.^ = L, D). However, they are below the experimental 
bounds even in the case of maximal mixing. 



In principle we could now proceed to impose bounds on neutrino masses [25 



4.35 eV, 
m Vn < 160 keV , 

77v ^ 23 MeV . (3.7) 



and lepton mixings from the partial Z width []T8 



r(Z^e ± /i T ) < 0.6 x 1(T 5 
T(Z^e ± r T ) < 1.3 x 1(T 6 



r(Z^iiT^) < l-9xl0" b , (3.8) 

in order to obtain constraints on the i?-parity breaking parameters^. This route has been 
taken by various other authors in models with spontaneous i?-parity breaking or in models 
with i?-parity breaking Yukawa couplings 0] and very recently also in the model under 



investigation here[[2S 



However, our goal is slightly more ambitious. Rather that trying to rule out a certain 
region of the seven dimensional parameter space we are more interested to see whether our 
predictive model with only three i?-parity violating parameters can provide the solution 
to actual problems. In particular, we want to find out whether the pattern of neutrino 
masses can provide a natural framework for a solution of the solar |]16| and atmospheric 



neutrino puzzles. The masses and mixing angles needed to solve the atmospheric neutrino 
puzzle are[]TT| 



ml r -ml ~ 1(T 2 eV 2 , sin 2 29^ ~ 1 . (3.9) 



For the solar neutrino puzzle there are three possible regions in parameter space |jl6|j,p9 

ml -ml ~5x 10~ 6 eV 2 , sin 2 
ml -ml ~ 5 x 1(T 6 eV 2 , sin 2 

.2 2 _ ln -10 „\ T 1 „:„2 



777 2 - 777 2 e ~ 10" iU eV 2 , sin 2 2^ M ~ 1 . (3.10) 



The first two regions correspond to an MSW solution [g(J and the third one corresponds 
to long- wavelength oscillations (LWO) [f2~9|l . Our analysis is similar to the one performed 
in ref. and |TB[ in models with spontaneous i?-parity breaking and in ref. ]3T[] for 



models with i?-parity breaking Yukawa couplings. The model under consideration here 
distinguishes itself from the later models by permitting only three i?-parity violating 

2 A much stronger constraint of ra VT < 100 eV can be obtained by requiring that the neutrino density 
does not overdose the universe [E7|. 
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m [in GeV] tan Q l 

Figure 3.6: Scatter-plot of 6 x 10 4 models with maximal mixing vs. (a) logio (mo/Ge V) 
and (b) tanfli for sin 2 2^ = 0.008 



parameters in addition to the usual minimal SUGRA parameters and is thus considerably 
more constrained. On the other hand, it does not predict any new particles below the Z 
mass and cannot easily be tested in present [28] or future collider experiments. 



The desired values of m 2 ^ — m 2 , sin 2 and sin 2 29 eVfi can be obtained by adjusting 
tan#i, tan #2 an d tan# 3 , respectively. Clearly the desired value of m Vr [eq. |3.9|| is some 
four to ten orders of magnitude smaller than the ones obtained for maximum .R-parity 
violation [fig. |3.5|| implying that tan#i has to be rescaled by some two to five orders of 
magnitude. Similarly we can fix sin 2 29 llVr and sin 2 29 ev by tuning tan #2 and tan #3. The 
value of m 2 ^ — m 2 c is then predicted as a function of the SUSY parameters. Note that the 
source of the lepton number violations and lepton flavor violations are soft mixing terms 
that are uncorrelated with the Yukawa couplings. Thus, there is no relation between 
the neutrino mixing angles and the CKM-matrix of the quark sector such as in SO(10) 
SUSY-GUT models and we should forget any prejudice towards small mixing angles. 

We begin our numerical analysis for the case of small mixing ( i.e. sin 2 29 eu = 0.008). 
In fig. p.6| we present the prediction of m 2 ^ — m 2 e obtained by scanning over the entire 
SUSY parameter space. The dashed line indicates the value required to solve the solar 
neutrino puzzle. The SUSY parameters are chosen as described in the last section. We 
present the prediction of m 2 ^ — m 2 e as a function of (a) m and (b) tan#i. The value 
corresponding to a MSW solution of the solar neutrino problem is indicated by a dashed 
line. We see that the value of m 2 — m 2 e is uncorrelated with mo and grows with tan#i. 
By requiring m 2 — m 2 = 5 x 10~ 6 eV 2 we find that 2 x 10~ 5 < tan#i < 5 x 10~ 3 . 



In fig. |37?] we present a histogram of the number of models iV that yield a particular 
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Figure 3.7: Histogram of the number of models that yield a particular prediction for 
m 2 — m 2 and a scatter plot of all the model in the tan (3— tan d\ plane. We use the same 



models as in fig |3.6| . The solid curve in (b) shows the upper limit on the i?-parity violating 



parameter from cosmology 

prediction for m 2 — m 2 e . We have used the same set of models as in fig. p.6| . The dashed 
(dotted) line indicates the value required to solve the solar neutrino puzzle via the MSW- 
effect (LWO). We see that the predicted value of m v in most models is about one order 
of magnitude below the value needed to solve the solar neutrino problem (the value of 
m Ve is negligible in all models). However, the number of models that predict the right 
value is still so large that no particular fine-tuning is needed. 

In fig. |3.71 (b) we have singled out all the models with 3 x 10~ 6 eV 2 < m 2 ^ — m 2 e < 7 x 
1CT 6 eV 2 . We find a clear correlation between tan (5 and tan 9± that allows to predict tan 9± 
to within a factor of three for a give value of tan (3 independent of the SUSY spectrum. Let 
us now compare the obtained value for tan^ with experimental constraints. The most 



solid constraints come from collider experiments. However, they are very weak||28| and do 
not constrain the region in parameter space we are interested in. Stronger constraints can 
be obtained by requiring that the baryon asymmetry not be washed out at the electro- 
weak phase transition 0. It means that at least one of the three lepton has to stay out of 
chemical equilibrium, ie. the lepton number violating Yukawa couplings for the electron 
have to satisfy 

* = =610 _ 4 tanfl lS mfl 2 sinfl3 < 1(J _ r (3 . n) 

V 2m w Ho cos p 

Note that eq. fj. 1 1| only constraints the product tan Q\ sin 62 oc itl Vt sin 62 oc m„ M . In 
fig. p. 7] (b) we see that the small mixing solution to the solar neutrino problem is compatible 
with eq. [3.11| (soid curve) in particular in the region of large tan j3. 
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Figure 3.8: Same as in fig |0| except sin 2 2d eVfi ~ 1 and (b) m 2 ^ — m 2 e = 510 5 eV 2 and 
(c) m 2 — m 2 e = 10" 10 eV 2 . Note that maximum v e —v^ mixing leads to much stronger 
constraints on tan 9\ 
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We will now investigate the region of large v^-v^ mixing i.e. we adjust tan 63 such 
that sin 2 29 eUfi ~ 1. This case yields a spectrum of squared mass differences very similar 
to the previous case. The result is summarized in fig. [3.8| . We see that the value of 
m 2 ^ — m 2 e needed for an LWO explanation of the solar neutrino puzzle lies very close to 
the maximum. It is obtained by roughly ten times as many models as the value needed 
for an MSW explanation. Of course it is impossible to draw any rigorous conclusion from 
such a scan since we have no objective way of weighing the probability for a particular 
set of parameters. However, in fig. |3.7| (b) we see that the large mixing MSW solution 
is incompatible with cosmological constraints eq. |3.11| (soid curve) while for some sets of 
parameters the large mixing LWO solution is allowed. 



3.3 LSND and hot dark matter 

In the last section we answered the question whether the solar and atmospheric neutrino 
problems could be solved simultaneously within the frame-work of our model. However, 
there are several other experimental hints that indicat the existence of neutrino oscil- 
lations. Maybe the most significant result arises from \x decay at rest which requires 



m 1 ~ m t e ~ 0-3 eV 2 and sin 2 29 eVii ~ 0.004|32[]. We perform a similar scan over the pa- 



rameter space as in fig. |3.7| and fig. |3.8|. Note that here the solar neutrino deficiency can 



be explaine by v e -v T oscillation. Thus, we fix tan#i by tuning m VT — m Ve = 5 x 10 eV 
tan 9 2 by tuning sin 2 29 ev ~ 0.004 and and tan 9 3 by tuning sin 2 29 eVr ~ 0.008. 

Our result is summarized in fig. |3ll The dotted curve shows the lower limit on m 2 — m 2 



compatible with the LSND result. The short dashed curve shows the upper limit on 
m v from collider experiments ||18|| while the long dashed curve shows the upper limit 



on the heavies neutrino mass from requiring that the neutrino relic density does not 
overdose the universe ^7|. We see that for large values of tan/5 there are a few sets of 
parameters that can simultaneously explain the solar neutrino deficiency and the LSND 
result. However, since in this scenario m v is quite large in comparison with m„ r in the 
case of the atmospheric neutrino puzzle it is clear that the i?-parity violating parameter 
is larger than its cosmological upper limit 0. 

One of the disadvantages of models with broken .R-parity is the instability of the light- 
est supersymmetric particle which cannot serve as a cold dark matter candidate. One 
the other hand, it has been suggested that the existence of hot dark matter (HDM) 
consisting of massive neutrinos with a sum of all the masses of about 7 eV would be 
desirable |33|]. Clearly, this condition is only compatible with both a solution to solar and 
atmospheric neutrino puzzle if all three neutrinos are almost mass-degenerate. This is 
incompatible with the hierarchical spectrum of our model (ie. m Ve <C m v ,m UT ). How- 
ever, it is clearly possible to solve the solar neutrino puzzle via v^-v^ oscillation while 
fixing m VT = 7 eV. E.g, take all the models from fig. |3.7| (b) and fig. |3.8| (b)-(c) and replace 
tan0a -> (7 eV/0.1 eV) 1 / 2 tan^ 1 and tan# 2 -> (7 eV/0.1 eV)" 1 / 2 tanfl 2 - This rescaling 
only changes the value of m Vr while the values of m v , m Ue and the i?-parity violating 



coupling constrained by eq. 3.11 remain the same. 
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Figure 3.9: Histogram of the number of models that yield a particular prediction for 
m v ~ m l e • We use the same set of parameters as in fig |3.5| except that we only consider 
tan/3 = 2 and tan/3 = 40. We fix tan 9 2 by tuning sin 2 26* ei , M ~ 0.004 and tan #3 by tuning 
sin 2 26 eUr ~ 0.008 
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Chapter 4 
Conclusions 



We have investigated the spectrum of neutrino masses and mixing angles in the MSSM 
with broken .R-parity. We have focused on the model where -R-parity breaking arises 
explicitly from dimension 2 terms in the superpotential. This model is characterized by 
only three additional parameters and can be embedded in SUSY GUT model without any 
constraints from the proton life-time. In this model, the -R-parity breaking is characterized 
by tan#i, the /i-t mixing by tan 6*2, and the e-/i mixing by tan# 3 . By adjusting these 
parameters we can solve the solar and atmospheric neutrino puzzles without fine-tuning. 
In particular, we have shown that a large hierarchy m VT ^> m v ^> m Ve is quite natural 
even in the case of maximum mixing and may favor a LWO explanation of the solar 
neutrino puzzle. To obtain the correct neutrino masses we need tan#i ~ 10~ x (x = 2 ~ 5; 
this ratio may be explained as the ratio of Mqu T /M p ) but no additional particles and 
no new intermediate scale. Requiring that the baryon asymmetry not be washed out at 
the electro-weak phase transition rules out the large mixing MSW solution, while some 
regions of the LWO solution and small mixing MSW solution are allowed. 

Finally, we demonstrate that there are regions in parameter space where the solution to 
the solar neutrino puzzle is compatible with either the LSND result or the existence of a 
7 eV HDM neutrino. 
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Appendix A 



The RGEs for the imparity violating 
terms 



Here we present the RGEs for the R symmetry breaking terms. They can easily be derived 
from ref. In our model without i?-parity violating Yukawa couplings and with diagonal 



lepton Yukawa couplings the RGEs for fij reduce toQ 

32tt 2 ^ = fj, (try L2 + N c try D2 + N c try U2 -3g 2 -g r 

32vr 2 ^ = w (yfayfa + N c try U2 - 3g 2 - g' 2 ) , (A.l) 

where we do not sum over indices in brackets. Furthermore, for our purposes it is more 
convenient to write Bj = Bjjfij with 

167r 2dSoo = ^ ALyL + ^ ADyD + ^ AUyU + 3Mg2 + M , gl2 ^ 



dt 

16vr 2 ^ = Af u) y^ )+ N c tTA u y u + 3Mg 2 + M'g' 2 , (A.2) 



1 After completion of this paper, a complete set of RGEs was presented in ref. j35| 
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Appendix B 

The mass at tree-level 



In our model the electro-weak symmetry is broken via a the vacuum expectation value 
(VEV) of the neutral CP-even Higgs boson fields and the sneutrino fields, v = (H) 
and Vi = (Li). We find it convenient to present all the masses and interactions in 
the basis where the interactions are -R-parity conserving and i?-parity breaking is only 
parameterized by the mass parameters /ij, B^, and after electroweak symmetry breaking 
also by Vi. Furthermore, it is convenient to define the unit vector uj = vj/v. 



In our notation we use calligraphic letters to denote the mass matrices and interaction 
matrices in the electro-weak basis and roman letters for the mass matrices and interaction 
matrices in the basis of the mass eigenstates. The hypercharges are listed in table IB? 



In models with broken .R-parity the neutrinos mix with the neutralinos. We can write the 
potential as 

V =±(ij )M x0 (ij ) T + R.c, (B.l) 
where ^ = i—iX', — iX 3 , ipjj, ^\ )[] and the 7x7 neutrino/neutralino mass matrix 



M x o 



( M' 



\ 





M 

symmetric 



m z s w Sf3 
-m z c w Sf3 




-■m z s w Ci3Ui \ 
m z c w ci3Ui 
-Hi 
/ 



(B.2) 



where we have defined c w = cos# w = m w /m z and tan/3 = v/v, s@ = sin/3, etc. We now 
define the unitary matrix Z such that the matrix m x o = Z*M x oZ~ 1 is diagonal and the 
value of the diagonal elements is ordered by magnitude. 



The mass matrix in eq. [B.2| still has two massless eigenstates corresponding to two neu- 



trinos. However, one neutrino acquires a mass 



m U3 ~ 



c} (Msj + M'cj) m\y? sin 2 ^ - e[) 
M'Mifi - 2s p cp (Ms* + M'cl 



(B.3) 



1 Here, the superscripts are the SU(2)l indices. Note that we have switched the third and fourth 
component as opposed to other authors fed for convenience. 
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Q 


U 


D 


L 




H 




i 

3 


4 
3 


i 

3 


-1 


-2 


1 



Table B.l: The hypercharges of the Higgs fields and the matter fields. Note that 
Y&c = -1% for $ = U,D,E 



at tree-level. Furthermore, the charged leptons mix with the charginos. In complete 
analogy to the MSSM we can write 

r = (-iX-iti.,) (B.4) 

The potential is given by 

V = fM x ±^ + + H.c. (B.5) 
The 5x5 mass matrix is given by 

/ M v^nvsin/? \ 

M x ± = V2m w cos (3u no -Vifi , (B.6) 

V V2m w cos (3ui fit yfjVo J 

We now define the unitary matrices U R , U L such that m x ± = U R * .M x ±(£/ L ) _1 is diagonal 
(in the MSSM U R = U and U L = V). Then we can write the chargino mass eigenstates 
as four component spinors 



x - ' 



Let us consider the Higgs sector. In the case of explicitly broken i?-parity (i.e. Hi 7^ 0) 
we find that the Higgs multiplets mix with the lepton multiplets. In particular, the 
sneutrinos acquire a non-zero VEV. The system of five equation corresponding to the 
minimum conditions of these fields can be solved most easily numerically by iteration. 
The tree-level mass matrix of the neutral, CP-even Higgs/slepton mass matrix is given 
by 



-Bf - m\spCfiUi hiVj + m L u + l m z[ c U 2uiU J + ~~ 

and tree-level mass matrix of the neutral, CP-odd Higgs/slepton mass matrix is given by 



M = ^? + < + X(4-4) - b j 

A V -B? HiVJ + rnl^ + lmlicl-sl^u, 

H.ml{ • ™0<*»P"A, (B.9) 

z V sin p cos puj cost pujUj / v ' 
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Finally, the two charged Higgs bosons mix with the charged sleptons. The mass matrix 
in the electro- weak eigenbasis if+ = (H + , Ef + , E R+ ) with a = 1, 8 becomes 

where = (L~)*, E R+ = E c where the submatrices are given by 

12 / 1 -sin2pX/ 
^ 2 w V-sin2/3M 7 2cos 2 /3m 7 m j -cos2/3(5 / j 

w V sin /? cos pu/ cos A puju, j J v y 



with Mjj given by 



/If 2/2 _ _ M y2 

m 00 — m j — m i0 



and 



m% = y^vl + y^vm, (B.12) 



M\ R = I AfjVi - Hiy^v I (B.13) 
-Af j v + fx y^v i 



M 2 RR = M 2 Eij - \Y E cmlsl cos 2(3 + y^y&l + VuVijVkVi (B.14) 

Note that the gauge dependent pieces arise from the gauge fixing terms in the gauge. 
They give masses to the goldstone-bosons m 2 A0 = ^ z ml and m 2 H ± = £w^w> respectively, 
Inside our one-loop calculation of Appendix D these unphysical fields are treaded just like 
the physical fields. 
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Appendix C 
Vertices 



In the following appendix we summarize a set of vertices in a notation that is convenient for 
SUSY theories with broken R parity [|28fl,pl|l. Note, that due to the minus signs in eq. 2.6 
the fields ip~ defined in eq. B.4 are minus the fermionic partners of H~ (x = 2,3,4,5) 
defined above eq. B.10. 

The Feynman rules are given by i r y l ^(V L L + V R R), -i(V L L + V R R), -(V L L + V R R), 
for couplings of a fermion pair to a gauge boson, a complex or CP-even scalar, and a 
CP-odd scalar, respectively (here, V L and V R is the generic notation for all the vertices 
listed below. Note that we have included the gauge coupling constant in the definitions 
oiV L ' R ). 



The W XnXt vertices are given by (n = 1, .., 7 and x — 1, .., 5) 

V2 



1 



®nx - 9 S n2 5 x i 7=S n3 5 x 2 , 



>na+2$xa ■ (C.l) 

The ZXxXy vertices are given by 



,L = -^diag(2si-2,2s^-l,2s^ (ii) 

O lR = ^-diag (2s^ - 2, (2s^ - , (C.2) 
no summation over indices in parenthesis. The Zx^Xn vertices are given by 

0" L = -i-diag (0,0,1, , 

0" R = -0" L (C.3) 
The H^XyXt vertices are given by 

1 ( gSi y 6 z2 for x = 1 

Qx yz = 7f x 9W,i + Vy-2 Z -2 for X = 2 (C.4) 
V { gS xy 5 zl - 2/JL 2z _ 2 <^2 for x = 3, 4, 5, 
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and Q* yz = Q xzy - The A° x x v xt vertices are given by 



-■xzy 



^xyz 



Qz,. z for x — 1 

xyz 
^xyz 



-QLz for x = 2, 3, 4, 5, 



and S* g = S% zy . The H°x° m Xn vertices are given by 

| [Smx+2 (8n2 ~ tan0 w 5 n i) + 8 nx+2 {8 m2 - tan 9 w 8 ml )] for X = 1 



Q' 



//P 

xmn 



f [<W+2 (<5„ 2 - tan 6> w 5„i) + <5 n:r+ 2 (<5 m2 - tan 6> w 5 m i)] for x = 2, 



where P — L,R. The vertices are given by 



^xmn 



S'xmn — | (#n2 ~ tan0 w 5 nl ) + 5 n:E+ 2 (5 m2 — tan# w 5 ml )] 



The if a XnX^ vertices are given by 



IL 

anx 



' YH^8 n i8 x2 + -^8 n2 8 x2 + g8 n3 8 xl for a = 1 

for a = 2 
for a = 3, 4, 5 
. for a = 6, 7, 8 



Vn-Ax-2 
~Ua-2x-2^n4 



Q 



anx 



'0 for a = 1 

-^L^nAa + ^ n2 £xa ~ 0<Wh2&cl f° r a = 2, ..,5 
. -yn-4a-5^2 + ^-2a-5 5 n4 for a = 6, 7, 8 



The u a UbXn vertex is 



~^abn 



The d a dbXn vertex is 



= Vab^n3 ~ -J=Yu9'8a-3b8nl 



Kbn = Vba-3 5 n3 + ^g8 ab 8 n2 + -^Y Q g'8 ab 8 nl 



V, 



dL 
abn 



"aim 



The u a df,x x vertex is 



The d a Ubxt vertex is 



— Vab^m — ~^Y D g' 8 a s b 8 nl 

= y ba -3$n± ~ -^g8 ab 8 n2 + -^=Y Q g'S ab S nl 



KbxC = -V D ab 8x2C 

VllC = (g8 ab 8 xl - yl_ 3 8 x2 ) C 



v ab x - -y a b°x2 

Kt = gS*b5xi - y ba - 3 8 X 2 
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Definition and more details on the charge conjugation operator C can be found in ref. 
The vertices in the mass eigenbasis are given by 



O p 


7 TT P (1 p 

^ nn^ xx^ nx 





ttP ttP ml P 

^ xx yy xy 





ry ry ff\ll P 

^mm^nh^fhn 


s 

xyz 


— // // // tS 

xx yy zz xyz 


OH 

xyz 


— J) J J J J 9 

xx yy zz xyz 


Qxyz 


= U H -°U R -U L ~Q^-~ 

^ xx ^yy^zz^xyz 


R 

^xyz 


= U H -°U L -U R Q R ~ 

xx yy < zz-^xyz 


o' p 

^anx 


- JJ H± 7 -JJ P O lP 

~ u aa n nn < ~> xx^afix 


o" p 

^5 xmn 


TI H ° 7 7 o" p 

u xx ^rnm^nn^xmh 


QllP 

^xmn 


TjA° y ry q/fP 


ydP 

' aax 


= u s ~u p v? p ~ 

^ aa^ xx r aax 


iruL/R 
'aax 


tju tjR/L^uL/R 
^ aa^ xx * aax 


V qP 

aam 


- TT^ 7 V qP 

u aa* J rnm * aam 



(q = u,d) 



(C.14) 



for P = L,R (here U L = V, and U R = U). Furthermore, we have assumed that the 
Yukawa matrices y u / D l L are diagonal. The unitary matrices U** diagonalizing the squark 
mass matrices are defined analogous to eq. 2.15. 
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Appendix D 

The One-Loop Corrected 
Neutrino /Neutralino 
Mass Matrix 



In this appendix we present the complete one-loop radiative corrections to the full neu- 
trino/neutralino mass matrix. Note that for the radiative generation of the neutrino 
masses only the diagrams involving Higgs fields are relevant. We have regularized the 



divergences by dimensional reduction |23|. The calculation has been done in the t'Hooft- 
Feynman gauge (£ = 1). Note that the self-energy diagrams and as a consequence also 



the running DR masses as defined in eq. |3.1| are gauge- dependent. There are seven 
types of diagrams contributing to the self-energy of the Neutrino/Neutralino: charged 
Higgs/chargino, CP-even Higgs/neutralino, CP-odd Higgs/neutralino, up quark/squark, 
down quark/squark, H /± /chargino, Z/neutralino. Their results are 



-.V 



I ( y R i v 

2 \ nm ' '-' 
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The -B and i?i functions are defined by 

[B (q 2 ,m\,m 2 2 ),q^Bi(q 2 ,ml,m 2 2 )} = -lQin 2 J 
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